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Abstract 

We study the dynamics of initial nucleation processes of photoinduced structural 
change of molecular crystals. In order to describe the nonadiabatic transition in each 
molecule, we employ a model of localized electrons coupled with a fully quantized phonon 
mode, and the time-dependent Schrodinger equation for the model is numerically solved. 

We found a minimal model to describe the nucleation induced by injection of an 
excited state of a single molecule in which multiple types of intermolecular interactions 
are required. In this model coherently driven molecular distortion plays an important role 
in the successive conversion of electronic states which leads to photoinduced cooperative 
phenomena. 

1 Introduction 

As the recent advance in laser technology has made it possible to generate arbitrarily designed 
optical pulses, control of quantum mechanical states of materials with those laser pulses has 
become to one of the central interests of the research of future device applications. For example, 
study on quantum information technology pQ has opened up a new aspect of device applications, 
i.e., utilization of the phase (coherence) of the quantum mechanical states. In those new fields, 
contrary to the conventional control methods of electronic/vibrational states, it is required to 
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keep the coherence of the internal states of materials during control processes, which means 
that we need to understand the dynamics of quantum-mechanical states in coherent regime. 

On the other hand, it was also found in various materials that injection of photoexcited 
states induces cooperative phenomena regarding with the change of structural, magnetic, or 
ferroelectric properties [2j[3], IU [6]. These photoinduced cooperative phenomena are considered 
to have a common mechanism, and many experimental and/or theoretical studies have been 
presented to make it clear [SJ Ej. in particular, when we are interested in controlling such 
cooperativity by designed optical pulses, it is necessary to understand the transient properties 
of the photoinduced cooperative phenomena more deeply. Thus, theoretical studies focused on 
the dynamics of the photoinduced cooperative phenomena are important at the current stage. 

As we have shown in the previous papers [TU| [TTj . nonadiabaticity of electronic transitions 
is a key to understand the temporal behavior of the above-mentioned phenomena. Dynam- 
ics of nonadiabatic transitions has been studied since the pioneering works by Landau [12J 
and Zener[13], and the bifurcation rate of wavefunction was analytically obtained in general 
cases[T3]. These studies mainly focused on the wavefunctions before/after nonadiabatic tran- 
sition, and hence the time evolution of wavefunctions itself is out of their scope. On the other 
hand, the dynamics of nonadiabatic processes has been considered to be important in, for ex- 
ample, photochemical reactions [T3], and hence computational methods of the dynamics have 
been proposed by many authors [16]. Since, however, those methods require the atomic coordi- 
nates/momenta to be treated as classical variables due to the limited computational capacity 
at the present time, they could discuss the wavefunctions after decoherence of the atomic de- 
grees of freedom takes place. As a result they mentioned the reaction yield or the absorption 
rate after various nonequilibrium processes. On the contrary, the initial nucleation processes 
in photoinduced cooperative phenomena involve consecutive switching of potential energy sur- 
faces (PESs) relevant to the dynamics of excited states, and the coherence of the wavefunction 
of electrons/atoms should be taken into account. In other words, the wavefunction at every 
moment should be pursued to understand the dynamics of the whole processes. Hence, not only 
the bifurcation rate of the wavefunction after nonadiabatic transition but also the wavefunction 
as a function of time is required, which means that the conventional methods are not suitable 
for the theory of photoinduced cooperativity. 

In this paper, we study the initial dynamics of the photoinduced domain growth, which is 
characterized by nucleation processes before decoherence of quantum-mechanical states takes 
place. We also focus on the photoinduced phenomena in which electron-phonon interactions 
play an important role and propose a minimal model to describe such cases. 

The organization of the paper is as follows: in Section [2] the molecular model is introduced 
and the method of calculation is described. In Section [3] the calculated results are shown. 
Section @] is devoted to discussion and conclusions. 
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2 models and method 



As we discussed in our previous papers[TOl [TT], nonadiabatic transitions between quantized 
states are particularly important to study the dynamical aspects of photoinduced cooperative 
phenomena. 

In this paper, we focus on the initial dynamics of a photoexcited state in interacting 
molecules, fully quantizing the relevant vibration modes. However, the dimension of the Hilbert 
space for the whole system increases drastically by quantizing atomic variables, which means 
that numerical calculation on those systems requires lots of computational resources. Thus, we 
employ a simplest model which is sufficient to describe the photoinduced nucleation processes. 
In the present model we consider molecules arrayed on a square lattice. Electrons relevant to 
the nonadibatic transitions are assumed to be localized in each molecule, and two electronic 
levels coupled with a single vibration mode is taken into account per molecule. The diabatic 
PESs with respect to the electronic states in each molecule cross with each other and that the 
nonadiabaticity in the dynamics is taken into account via "spin-flip" interaction between two 
electronic states. This model is known as a simplest model to discuss the relaxation dynamics 
of, e.g., photoisomerization of molecules [T7J. As for the intermolecular interaction, we take into 
account vibrational coupling and the Coulomb interaction between excited state electrons. j3 
affects to induce molecular distortion by the excited electrons in the adjacent molecules which 
is also of the same order as the other interaction terms. 

Hence, the Hamiltonian in the present study is described by: 

y + — ^ + (V2hLU 3 sq? + efiuj + s 2 hco)n^+ \a r x \ 
+ [a^ 2 (uf - (3n?)(up - /3hp) - {V - W(u? + up)}h?np\, (1) 

where and up are the momentum and coordinate operators for the vibration mode of a 
molecule at site r, respectively. The second sum which gives the intermolecular interaction 
is taken over all the pairs on nearest neighbor sites, where the Coulomb interaction between 
excited state electrons are modified by molecular distortion. The vibrational period of an 
individual molecule is denoted by T = 2n/u in the rest of the paper. 

A schematic view of the present model is shown in Fig. [TJ The two diabatic PESs for 
an individual molecule are crossed with each other, and the nonadiabatic coupling constant 
A acts to separate them into two adiabatic PESs. We chose the values of the parameters as: 
£ = 1.6,5 = 1.4,1/ = 1.1, W = 0.2, a = 0.1, p = 0.2, and A = 0.2. Although those values 
are typical for organic molecules as for electron-vibration coupling[18] and the intermolecular 
Coulomb interaction [17], the other parameters are not easy to determine their values either from 
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theoretical calculations or experimental results. We only mention that the order of magnitude 
for the parameters is estimated referring to those for typical organic materials. The electronic 



[ V I I I I 

w 

a 
P 

Figure 1: Schematic view of the model. Circles denote the molecules with two electronic states 
and a vibrational mode. Adiabatic potential energy surfaces for an individual molecule is shown 
in the inset. 

states at site r are denoted by | |)f (ground state) and | |)f (excited state) and a\ [i = x,y,z) 
are the Pauli matrices which act only on the electronic states of the molecule at site r. 
denotes the density of the electron in | \ )? which is rewritten as n r = a r z + 1/2. The model and 
the notations of the parameters are schematically shown in Fig. [TJ 

Details of the quantization procedure of the vibration mode of each molecule will be pub- 
lished elsewhere. We mention that the basis set for the vibronic states is composed of the Fock 
states shown in Ref. [19]. The phonon dispersion relation of the vibration mode is given by 

f2(/c) = oo^Jl + 2a(cosk x + cosk y ), (2) 

where (k x , k y ) denote the reciprocal lattice vector of the square lattice, and the lattice constant 
is taken to be unity. The quantized states on each diabatic PES of a single molecule are the 
vibronic states \na)? (n = 0, 1, 2, a =], {) in the Fock representation, where the coordinate 
of the molecule is labelled by r . \n |) is related with \n |) by 

|nT) = |T)U|e s(ot+a) |n|), (3) 

where e s(at+a ) denotes the translation operator in the vibration coordinate space[20]. We note 
that this Ising-like model is similar to the one to study the thermodynamical properties of the 
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J aim- Teller effect [2T| . though the nonequilibrium dynamics of the excited states in the model 
has not been understood. 

We obtain the numerical solution of the time-dependent Schrodinger equation for the Hamil- 
tonian (TT|) by the Runge-Kutta method. In each series of calculations, one of the molecules 
on 128x128 lattice is initially in the Franck-Condon state, while the others are in the ground 
state, which corresponds to the injection of a photoexcited state to a single molecule at the 
origin. In solving the Schrodinger equation, we applied a mean-field approximation in which 
the contribution of the wavefunction at the nearest neighbor sites is substituted by the average 
value with respect to the wavefunction The detail of the approximation is described in 

the appendix. We only mention here that this approximation is equivalent to decomposing the 
wavefunction of the total system |<I>(£)) into a product of the wavefunctions at each molecule, 
i.e., 

\$(t)) = |0(t)) rl ® \<j>(t)U ® ."I0(*))r&, (4) 

where N denotes the number of molecules in the system. Thus, we can solve the differential 
equation for each molecule when only the average values of the properties for adjacent molecules 
are provided, which means that the present calculation method is suitable for parallel comput- 
ing. Hence, we have made it possible to handle more than 10000 molecules by the present 
method. 



3 calculated results 

We have pointed out that the population of the excited electronic state | |)r is suitable for 
understanding the dynamics of the initial nucleation processes [TOl ITT] . Thus, we show in Fig. 
H] the population of | t)f for 48x48 sites around the injected excited state ("seed") defined by 

N(r,t) = ($(t)\nrmt)), (5) 

for t — 0, 5T, and 10T. Figure [2] shows that the number of molecules in the excited electronic 
state increase surrounding the initially excited molecule. Those molecules will constitute a pho- 
toinduced domain observed in many experiments [2J El HI El E] and thus the present calculation 
described the initial processes of the photoinduced cooperative phenomena, %. e. photoinduced 
nucleation triggered by an injected excited state. 

Since the population transfer is induced by adiabatic transition on the PESs in Fig. [lj 
molecular distortion is also relevant to the nucleation processes. Hence, molecular distortion 
defined by 

c(r,t) = mt)\ u ?mt)), (6) 



5 




6 



is calculated to discuss different aspects of the nucleation processes. Figure [3] shows ((r,t) for 
t = 0, 5T, and 10T. We found that the molecular distortion is coherently driven in the system 
and that the vibrational energy is propagated by coherent phonons. 

The elementary processes of the nucleation is understood by comparing N(r,t) and C(r,t) 
as functions of time. Although these properties seem to have similar temporal behavior at first 
glance, we found that the difference in their initial growth dynamics is a key to understand the 
details of the nucleation processes. We show in Figs. E]-(b)-(f) N(f,t) and C(^) as functions 
of time, where f corresponds to the molecular sites shown by filled circles in Fig. IU-(a). As we 
mentioned, the boundary between the converted (excited state) domain and the ground state 
domain moves to extend the former in the system. Thus the values of N(f, t) and £(f, t) rise at 
later time as the distance to the "seed" increases. Comparing Figs. HJ-(b)-(f), we found that the 
distance to the "seed" should refer to the Manhattan distance in discussing the time of rising 
of those values when \ f\ is small. However, as the growth process proceeds, the number of the 
excited molecules increases and the conversion of the electronic states takes place as in the 
continuous systems and the boundary between domains becomes a circle (Fig. El-(c)). Hence, 
the Euclidean distance becomes appropriate for the distance between two molecules. This 
behavior is reminiscent of the propagation of coherent phonons, and the nucleation processes 
are driven by those coherent motion of the molecules. 

Figures H]- (b)-(d) also show that an oscillating component in both N(f, t) and £(r, t) appears 
before transition to | f)r is realized. We found that, after N(f,t) and/or ((r,t) begins to 
grow, it takes longer time to complete electronic state conversion in the molecules close to the 
"seed" than in the others. To be more precise, a precursor to the electronic state conversion is 
observed in molecules at the nearest neighbor and the next nearest neighbor of the "seed" . This 
behavior is understood as a preliminary process required to overcome potential energy barrier 
at the beginning of nucleation. Thus, in the present model, the first process of the nucleation 
corresponds to making a cluster of the ~ 10 converted molecules around the "seed." Then the 
other molecules will suffer electronic state conversion which smoothly takes place as Figs. HJ-(e) 
and (f) show. 

The above properties are relevant to the mechanism of the photoinduced cooperative phe- 
nomena, i.e., each molecule in the ground state receives energy through the vibrational coupling 
a first and the molecule begins to vibrate. Then, the nonadiabatic coupling A induces popu- 
lation transfer with the assistance of the Coulomb interaction V, and thus the molecule turns 
to belong to the converted domain. Once each molecule belongs to the converted domain, mo- 
tion of N(r,t) and C(^) almost disappear and it does not return to the ground state during 
the current simulation, since intermolecular interactions make it remain in the excited state 
domain. We also point out that the intermolecular Coulomb interaction V enhances the pop- 
ulation transfer rate particularly when the number of adjacent molecules in the excited state 
increases. 
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Figure 3: Molecular distortion ((r,t) on 48 x 48 lattice for (a) t — 0, (b) t — 5T, and 
t = 10T. 
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Figure 4: Time-dependence of N(f,t) and C(^)- for each figure corresponds to the filled 
circles in the top figure, and the "seed" is denoted by the shaded circle. 



9 



150 




5 10 

time(T) 



Figure 5: iV tota i(i) as a function of time. The dotted line proportional to t 2 ' 3 is drawn as a 
guide for the eyes. 



Figure H]-(b) also shows that the molecular distortion does not grow for ~ 3T after photoex- 
citation. Hence, it takes ~ 4T for the population transfer to be completed even for the molecules 
in the nearest neighbor sites. Since the population transfer reflects on the electronic/optical 
properties of the molecules, the change of those physical properties takes place ~ 400 — 800 
fsec (T is ~ 100 — 200fsec in typical organic molecular systems) after photoexcitation. 

The sum of the excited state population N tota \(t) = Y^r N(r, t) indicates the measure for 
the growth rate of the photoinduced domain. Figure [5] shows that iVt ota i(i) increases as ~ t 2 ' 3 . 
except in the very first stage of the domain growth. As a result we obtain that the radius 
of the photoinduced domain behaves as ~ t , which is understood by the picture that the 
growth of the domain is predominantly driven by propagation of coherent phonons rather 
than diffusion processes. As the vibrational coherence is lost, diffusion process becomes more 
important and the domain growth will slow down to make the radius of the domain increase as 
~ \fi. Since vibrational coherence survives for a few psec in typical organic molecules [18], the 
present calculation is valid only in the time range studied in this paper, and the decoherence 
of the vibrational states should be taken into account to study the growth dynamics of the 
photoinduced domain in a longer time scale. 
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4 discussion and conclusions 



In this paper we study the coherent dynamics of photoinduced nucleation processes in organic 
molecular systems. When a single molecule excited to the Franck-Condon state, it induces 
distortion of adjacent molecules, and the excitation energy is transferred to the other molecules 
coherently. Once the molecules start to vibrate, the electronic state conversion from | I)? 
to | f}^ takes place and thus photoinduced domain grows. This is the basic scenario of the 
initial photoinduced nucleation processes where coherent phonons play an important role. In 
fact, the size of the converted domain (diameter) is almost linearly increases as the nucleation 
proceeds, which shows that energy diffusion is subsidiary in the initial processes. However, as 
the decoherence of vibrational states takes place, excitation energy propagation in the system 
will be dominated by diffusion processes, and hence the growth rate will be oc \/i after all. We 
mention that these properties will be reflected on the time- resolved spectra of e.g., reflectance, 
absorbance, or Raman scattering intensity and that the ultrafast spectroscopy will give a key 
to understand the coherent nature of the nucleation processes. 

We stress that, contrary to the one-dimensional case [5], the domain growth by successive 
conversion of the molecules does not take place in higher dimensional cases if intermolecular 
interactions other than vibrational coupling a is neglected. This is a particular property in 
higher dimensional systems, since it was pointed out that the converted domain grows only 
by the intermolecular interaction between molecular distortion [8]. In general, the value of a 
should be less 1/M where M is the coordination number of the lattice. Hence, the maximum 
value of a is smaller as the dimensionality of the system is higher. The present study shows 
that, even in two-dimensional systems, a is not sufficient to induce cooperative phenomena, and 
thus other interactions such as Coulomb interaction between electrons are necessary. Since M 
is larger in higher dimensional systems, this result is applicable to three-dimensional systems, 
and thus we conclude that the photoinduced nucleation processes are realized by cooperation 
of various types of intermolecular interactions in general cases. We mention that the results 
described in this paper are independent of the model and the values of the parameters for which 
the nucleation takes place. 

In the present paper, we assume that only a single relevant vibration mode exists in each 
molecule. However, the nonadiabatic transition within a single molecule is strongly affected by 
the structure of the PESs. In particular, when multiple vibration modes are taken into account, 
the dynamics of the nucleation processes depends on the topological structure of the intersec- 
tions of the PESs, e.g., existence of conical intersections. Hence, ab initio electronic-structure 
calculations of specific materials are important for a detailed discussion of such material- 
dependent features of the nucleation processes, and the dynamics calculation in the present 
paper should be combined with those electronic-structure calculations in the future. We, how- 
ever, stress that the present results give the basic properties of the nucleation dynamics in 
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coherent regime and that the qualitative feature of the domain growth is sufficiently discussed 
in this paper. 

In order to control the nucleation dynamics by outer field, e.g., laser pulses, we should 
estimate the effect of decoherence of the quantum-mechanical states. We point out that it is 
possible to take into account the decoherence by embedding the system in a large 'reservoir' and 
by tracing out the dynamical variables regarding with the reservoir. We also stress that such 
studies will contribute to realizing the coherent control of the photoinduced domains, which 
will be important both from a physical point of view and device applications [22]. 
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A equations of motion for numerical calculations 

The Schrodinger equation for the system is written by 

i^|*(t)>=«|*(t)>. (7) 

When the wavefunction | <!>(£)) is decomposed by those for individual molecules as shown in Eq. 
(H]), each component \(j)(t))p obeys the following equation: 

d I t& uP'iA -\ 

ifo\<t>{t))? = | 2 ^ 2^ + W^ uj ^ S( 1t + + s 2 hu)h i ?+ \a r x > \<j){t))p 

+ J2l auj2 ( u r ~ Pn?)(up - Pnp) -{V - W(u? + up)}n?hp]\(j)(t)}?, (8) 
P 

where the sum is taken over the nearest neighbor sites of r. 

We apply a mean-field approximation to Eq. (jSJ), i.e., the effects of molecules in the nearest 
neighbor sites are substituted by their average values. Hence, after canonical quantization of 
the pjf and up, the following equations of motion for each molecule are derived: 

d 

i-Q^\4>{t))r = hu)[ala?+ {s(at+ a?) + e + s 2 }h?+ \a r x ]\(f){t)) ? 

+ [ auj2 b (4 + a f) ~ Pn?} ( {up )-P(np)) 
p 

- {V - WM4+ a?) + (up)}}h r {n p )] |0(t)>* (9) 



12 



where a}} denotes the annihilation (creation) operator of the vibration mode and 7 = sJU/2uj. 
Intermolecular interaction is taken into account through the average value over the decomposed 
wavefunction and 

{X P )?= MWlXMit))?. (10) 
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